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Introduction 
Aquifers, oil, and gas reservoirs, as a rule, have a heterogeneous structure at the micro- and macroscale, 
[1]. Heterogeneous reservoirs on a macro scale consist of different zones with different, sometimes very 
strong, filtration-capacitive properties, i.e. porosity, permeability, etc. Zones with well porosity and 
permeability are well conductors for liquids and various substances suspended or dissolved in fluids. A 
typical example of heterogeneous formations is fractured porous media (FPM), [2,3], the structure of which 
is represented as a system of fractures surrounded by porous blocks. 
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In [6] a transport model in a medium with double porosity was considered taking into account the 
reversible and irreversible deposition of colloid particles in both zones and the first-order equilibrium mass 
exchange between the zones. In each zona, i.e. in fractures and porous blocks, a reversible and irreversible 
deposition of particles with various characteristics occurs, described by linear equations. An analytical 
solution to the problem is obtained, which is used to describe the results of previous experiments, [5]. 
Coefficients of mathematical models are defined as the solution to the coefficient inverse problems (CIP), 
known as identification problems, [8]. It is assumed that the coefficients of the equation depend on the 
spatial coordinates and are independent of time. The statements of the problems are based on the use of 
uniqueness theorems for the solution of the CIP proved in [7], [11], [12], [15]. To obtain a unique solution 
of the CIP, it is required to set an overdetermined set of boundary conditions on the boundary of the zone: 
the function for which the equation is written or its normal derivative. 

Coefficient inverse problems (identification problems) have become the subject of intensive study, 
especially in recent years. Interest in them is caused primarily by their important applications. They find 
applications in solving problems of designing oil reservoir development (determining the filtration 
parameters of reservoirs), [10,10,14,16,17,18] in solving problems of environmental monitoring, etc. The 
standard CIP statement contains a residual function, which depends on the solution of the corresponding 
problem of mathematical physics, [14]. Methods for numerical solution of CIP in connection with their 
applications in underground hydrodynamics were developed in [7], [8], [9], [11], [13]. 

In this paper, an inhomogeneous two-zone medium is considered a single-zone medium with a source 
(sink). The second zone is modeled through the source (sink). This approach is fundamentally new because, 
in fact, the bicontinual medium is presented as monocontinual one. The validity of this approach is justified 
by the convergence of the results on the basis of the mono-continuous approach to the corresponding results 
of the bicontinuous approach. In the work, this is done by minimizing the residual function. In addition, it 
is assumed that in both parts of the first zone there is reversible adsorption of particles with the 
corresponding kinetic equations. Identification of parameters in the source (sink) term in the mass balance 
equation is carried out by solving the corresponding CIP using data from [4]. 


First zone 


Fig. 1. Scheme of solute transport in a two-zone medium 


2. The mathematical model and its numerical implementation 
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An inhomogeneous porous medium is considered, consisting of well-permeable and relatively low- 
permeable zones, the diagram of which is shown in Fig. 1. The parameters in the first zone are indicated by 
index 1. There are two sections in zone 1, in each of which the particle deposition with reversible 
nonequilibrium nonlinear kinetics occurs. It is believed that such processes also occur in the second zone, 
but we will not write equations and conditions for it. With the second zone, there is an exchange of 
substance, which we model by the fractional-order derivative in time of the solute concentration in the first 
zone. Consequently, in contrast to [4], the concentration field in the second zone is not considered. Note, 
that the fractional approach was previously used in [19], [20], [21]. 

The equations of solute transport in one-dimensional case are written as 
as, oS, ,c ac arc, eC, 
p—+p + 0, +a, —Ov,—., (1) 

Ot Ot Ot ” Ox 


where f is time, s, x is distance, m, D, is longitudinal dispersion coefficient, m i, S,v, is the fluid 


= 0, D, 


Ox? 


velocity, m/s, C, is volume concentration of the solute in the fluid, 5. and ee are concentrations of 
deposited particles, m*>/kg, @ is porosity, m*/m’*, ( is medium density, kg /m°*,a, is retardation 
factor related to the mass exchange between two zones, 5’ es y is the order time derivative with respect 
totime, O0<y <1. 


The deposition of particles in each of the sections of the first zone is reversible with the difference 
kinetic equations 


as. 

P a =6K,,C, FIP, as (2) 
os 

P= OK ON haa a2 (3) 


where k_,, K_, are coefficients of solute deposition from the fluid phase to the solid phase, s', K.,,,; 
k are coefficients of substance detachment from the solid phase and transition into liquid, 5’. 
ad2 p q 


Let a fluid with a constant solute concentration C, be pumped into the medium initially saturated with 


pure (without particles) liquid from the initial moment of time. Let us consider such time periods where the 
concentration field does not reach the right boundary of the medium, x = oo. Under the noted assumptions, 
the initial and boundary conditions for the problem have the form 


C,(0,x)=0, S,,(0,x) =0, S,,(0,x) =0, (4) 
C,(4,0) =c,, (5) 
es (t,00) =0. (6) 
Ox 


The problem (1) - (6) although linear, obtaining an analytical solution is difficult, because three 
concentration fields must be found at the same time. Therefore, to solve the problem, we use the finite 
difference method. In the considered region (2 = {(t, x)O<t<T, O<Sx< oo} a uniform grid was 


introduced 


= T a —— 
On lt, ast, =], x; et Jt, J=9, it. 
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where / is a sufficiently large integer chosen so that segment [0, x], x,=ih, overlaps the area of the 


calculated change in the fields C7, Sai, and Sa2. h is the grid step in the x direction. 
In the open grid area 


T _ _ 
O4={ 8), = x, =th, a aa J, i=1,/—-1, 


equations (1), (2), (3) were approximated as follows 


ju j iH j JH j 
p Sal" Sy, Sad" = Gall yg GAG, 


x 
_ a T T(2-y) 
HC) =(C)' gen, CM -©)De" 
«| 5 a ( i, ((j k +1) ”—(j-k) are J : yr |- (7) 
= 0D, (C, il 2(C,); a (C, i+1 Oy, (CJ; 7 (Ci, 
h? h 
(Sa. ~ (S,,) j J+ 
P a = BEC): = PK a (Sai); ’ (8) 
(S,,)/" ~ (Ses) J jtl 
P 7 - AK. (C,); — PK 42 (S,,); ? (9) 


where (C,)!, eae (S..)? are grid values of functions C,(t, x), S,,(t, x), S,.(t, x) ata given 
point - a) 


From the explicit grid equations (8), (9) we determine (S ai er (8 a iy ‘. 


(SRS= Py(S)-+ Pros (10) 
(So) NGOa): + Das (11) 
where 
1 TOK. 
— ane a C te 
- Tit TK a1 wa P mle PTR 4, vi 
1 TO. k 
Gin ~ — (Ci 


The grid equations (7) are reduced to the form 


AC) =B(C TEC), = Fis (12) 
where 
Ae 0 Dt Ze 

h? h 


20Dt Ovt at’ 
,=6,+— ++ + 
h h TQ-/” 
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alle 
h? 
. a jt a 
(Fy) =(,+ re repo ee Cal al, Panag 


ip (Gaia =G209G), = G=erl)*sG-0 GC); 


The following procedure of computing is used. From (10), (11) (S ae (S Ps are determined, 
then we solve the system of linear equations (12) by Thomas’ algorithm in order to calculate (C )i 7“, 


Since P,.> 4, <1, schemes (10), (11) are stable, and for (12) the stability conditions of the Thomas’ 


algorithm are satisfied. 
To assess the performance of the proposed model, it is important to compare the results with the 
corresponding results, [4]. To do this, we compare the source (stock) terms a(C,—C,) in [4] and 


Y 


a - in (1). To quantify the proximity of the results based on the curves was calculated 


2 


6,=((,-1,) dx (13) 


for a given value of t, where L is the conditional boundary of the medium to which the concentration 
profiles extend, 
0’C 
L=af€;—Cs, 1, =a, —. 
Ot” 


The proximity of the terms 7, and J, should guarantee the proximity of the concentration fields C, 


determined using the proposed approach and the model, [4]. To estimate their proximity, we use the 
standard deviation (13), only for C, determined on the basis of two models, i.e. 


6, =(c ce} ax, 
0 


where Cre concentration field C, (t,x) for a given t, determined according to [4], and C © the 


same as defined here. 
For other moments t and @, d,, y different estimates can be obtained for 6, and 0,. In principle, 


to approximate the two models, it is important to set and solve the corresponding coefficient inverse 
problems by determining of a,, v for a given value of @ or, conversely, determining of a fora given a 


and 7. 


3.Numerical results and their analysis 


In the calculations following initial values of parameters are used: 
c,=0, v,=10*m/s, D =v,-a@,, @, =0,005m, p=1800kg/m’, 0, =0,1, 
k=3:10%s", ky =2,5:10°%s", — k,, =4-104s7, k,,, =2-10%s" and various a), ¥. 
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We minimize the functional 
TL 

O(a,y)=[[(j-Lyade (14) 
00 


that characterizes the standard deviation of J, from /, for the entire time period. The calculations 
show that the minimum value of ®(a,,) is achieved at a, = 0,0006 y =0,8. 

The proximity of the terms 7, and /, should guarantee the proximity of the concentration fields C,, 
determined using the proposed approach and model, [4]. For this, the corresponding profiles are plotted 
for the data, obtained through minimization of ®(a,,y) (Fig. 2). As can be seen from the graphs, the 
solutions are close to each other. 

For a numerical estimation of their proximity, we use the standard deviation of the type (14), only for 


the one determined on the basis of two models, i.e. 
TG 


F(a,,y)= [f(CO-C?F deat, 


00 
where Cc is the concentration field C,(t,x) for a given t, determined according to [4]. Cc is the 


same defined here. For the cases analyzed above, the following minimum value of F(a,,y) was obtained 
0,00234738654452 for a,= 0,0006 y =0,8. 


Fig. 2. Comparison of concentration profiles obtained on the basis of two models a7 =0,0006, 


Cig 


y=08, 


The analysis shows that the simpler model proposed here, with an appropriate choice of parameters, 
can satisfactorily describe the results of a more complex model, [4]. 
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